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We have investigated the temperature dependence of the H \\ c flux line lattice structural phase 
transition from square to hexagonal symmetry, in the tetragonal superconductor buNi2B2C (T c = 
16.6 K). At temperatures below 10 K the transition onset field, HiiT), is only weakly temperature 
dependent. Above 10 K, H2(T) rises sharply, bending away from the upper critical field. This 
contradicts theoretical predictions of H2(T) merging with the upper critical field, and suggests that 
just below the H C 2(T)-cwcve the flux line lattice might be hexagonal. 



Studies of the topology of the magnetic flux line lat- 
tice (FLL) in type-II superconductors have a long history. 
Early neutron scattering experiments on 1ow-k supercon- 
ductors such as niobium revealed a multitude of different 
FLL symmetries and orientations, mainly determined by 
the symmetry of the atomic crystal structure in the plane 
perpendicular to the applied field [Q . This is not surpris- 
ing since deviations from the hexagonal FLL character- 
istic of an isotropic superconductor, and the locking to 
the crystalline lattice are driven by the symmetry of the 
screening current plane and by nonlocal flux line interac- 
tions within a range determined by the coherence length 
£o, and the electronic mean free path £. Later, similar 
effects were also observed in the strong type-II supercon- 
ductor VaSi with Kwl7fl, demonstrating that nonlocal 
effects can be equally important in high-K materials. 

Over the last couple of years, effects of nonlocality 
have been clearly observed in the borocarbide supercon- 
ductors. The borocarbides are quaternary intermetallics 
with stochiometry ift^P^C (R = Y, Gd-Lu) and a 
tetragonal unit cell (14/ mmm) ||. These materials have 
attracted attention due to the coexistence of supercon- 
ductivity (R — Y, Dy-Tm, Lu) and antifcrromagnetic 
ordering (R = Gd-Tm). The borocarbides are strong 
type-II superconductors with Ginzburg-Landau (GL) pa- 
rameter, k — 6-15. The discovery of a square FLL in 
most of the H \\ c phase diagram which undergoes a 
smooth transformation into hexagonal symmetry at fields 
below 1 kOe || , was the first observation of a purely field 
driven FLL symmetry transition. 

Using nonlocal corrections to the London model and 
incorporating the symmetry of the screening current 
plane obtained from band structure calculations, one is 
able to calculate the FLL free energy, and thereby to de- 
termine the stable FLL configuration in different fields in 
the borocarbides . The model succesfully describes the 



nature of the FLL square to hexagonal symmetry evolu- 
tion in the borocarbides with H \\ c as the applied field is 
reduced. Qualitatively, this can be understood as driven 
by the four fold basal plane anisotropy which makes the 
vortex current paths "squarish" close to the core. At 
high densities this leads to a square FLL, whereas at low 
fields the system appears isotropic resulting in a hexag- 
onal FLL 0. The onset of the transition occurs as the 
field decreases, commensing at a critical field H 2 , deter- 
mined by the range of the nonlocal interactions [^| . Later 
studies of the FLL in a set of samples doped with vari- 
ous amounts of impurities to systematically change the 
nonlocality range, showed a quantitative agreement be- 
tween the predicted and measured onset of the transition 
||[)| . The model also predicted a reorientation transition 
of the rhombic FLL at a field H\ <§C H 2 , which has been 
subsequently observed Jl0| , pd| . 

Whereas there is a growing experimental and theoret- 
ical understanding of the FLL geometry and transitions 
for fields well into the mixed state, there has been little 
attempt to understand what happens at fields just be- 
low the upper critical field. In this letter we report on 
detailed measurements of the temperature dependence 
of the structural transition, H2(T), in LUNI2B2C with 
the applied field parallel to the crystalline c-axis, over a 
broad range of temperatures and applied fields. Below 
10 K, H 2 (T) is essentially constant. Above 10 K, H 2 (T) 
rises sharply, eventually bending away from the H C 2(T) 
curve, suggesting the existence of a hexagonal FLL in the 
region just below the upper critical field. 

The FLL in LuNi 2 B 2 C (T c = 16.6 K) was studied us- 
ing the small angle neutron scattering spectrometer on 
the cold neutron beam line at the Ris0 National Labo- 
ratory DR3 research reactor. The sample used was a 1 
gram single crystal, grown from a high temperature flux 
Q using isotopically enriched n B to enhance the neutron 
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transmission. Incident neutrons with wavelengths in the 
range A n = 6 — 15.6 A and a bandwidth AA n /A n = 24% 
were used, and the FLL diffraction pattern was collected 
by a position sensitive detector at the end of a 6 m evac- 
uated tank. Magnetic fields in the range 1.5 to 10 kOe 
were applied parallel to the crystalline c-axis and to the 
incoming neutrons. In order to keep the resolution of 
the spectrometer unchanged, the neutron wavelength was 
varied in concert with the applied field to keep the scat- 
tering angle, 29 = gioA n /27r, constant. Here q w is the 
FLL scattering vector shown in figure |. The measure- 
ments were performed at temperatures between 2 and 20 
K, following a constant field cooling procedure. A back- 
ground signal obtained above T c was subtracted from the 
data. 

In figure |l| we show three different FLL diffraction pat- 
terns obtained in LUM2B2C. In quasi-2D structures such 
as that of the FLL, the real space flux line configura- 
tion and the reciprocal space diffraction pattern are re- 
lated by a simple 90° rotation and relabeling of the axes. 
The diffraction pattern therefore directly shows the FLL 
symmetry and orientation, with the caveat that scatter- 
ing from differently oriented domains is superposed. For 
the square FLL there is no orientational degeneracy with 
respect to the atomic lattice, due to the tetragonal crys- 
tal structure of the borocarbides, as reflected in the top 
panel of figure [l]. Here both the FLL (1,0)- and (1,1)- 
reflcctions are clearly visible, and resolution limited both 
in the radial and the azimuthal direction. The square 
FLL is aligned with the scattering vector qio along the 
[110] crystal direction. 

However, as soon as the FLL is distorted away from a 
square symmetry there is a two fold degeneracy, since ori- 
entations with the long rhombic unit cell diagonal along 
either the a- or the 6-axis have the same energy. Exper- 
imentally this will first manifest itself as an azimuthal 
broadening and, when the distortion increases as a split- 
ting of the (l,0)-reflections, with the two peaks each be- 
longing to a different domain orientation |5j . The magni- 
tude of the splitting will eventually reach 30° if the FLL 
is transformed to a perfect hexagonal symmetry. The 
azimuthal broadening is evident in the middle panel of 
figure 0, and in the bottom panel the (l,0)-reflections are 
clearly split. The effect of the transition is also reflected 
in the square FLL (l,l)-reflections. Ideally a radial split- 
ting should be observed, with one peak moving towards 
shorter q and gaining intensity and the other moving to 
longer q and fading, to become respectively hexagonal 
(1,0)- and (l,l)-reflections belonging to different domain 
orientations. However, due to the very rapid decrease 
in the scattered intensity with increasing q, one only ob- 
serves a radial broadening and a relative increase in the 
intensity. 

In general, the azimuthal splitting was only directly 
observable when above ~ 10°, comparable to experimen- 
tal resolution. To determine the transition onset field 



H2(T), we analyzed the azimuthal intensity distribution 
in each diffraction pattern at a scattering vector q = gio, 
where qw — ax 2n(B /(fro) 1 ^ 2 > 4>o — hc/2e is the flux 
quantum and a is a constant which depends on the FLL 
symmetry: With a FLL apex angle [3 as shown in figure 
||, a = (sin/3) -1 / 2 , and hence a = 1 for the square case 
((3 = 90°) and (2/V3) 1 / 2 = 1.075 for the hexagonal FLL 
(/3 = 60° or 120°). In the top panel of figure |^ we show 
the temperature dependence of the azimuthal width for 
H = 3 kOe, obtained by fitting to a single Gaussian. A 
clear broadening is seen for T > 10 K. The low-T width 
of 13.8° ±0.2° is a direct measure of the azimuthal resolu- 
tion, which was kept constant througout the experiment. 
To determine the azimuthal splitting, the intensity dis- 
tribution at each temperature and field was then fitted 
to two Gaussians with a width fixed of 13.8°, as demon- 
strated in the top panel inset. In the bottom panel of 
figure H the azimuthal splitting is plotted versus tem- 
perature for a number of fields between 1.5 and 6 kOe. 
At low temperatures the splitting goes to 0° ± 2° as the 
field is increased. As T is raised the splitting increases 
smoothly for H < 2 kOe, and abruptly above a certain 
threshold for H > 3 kOe. The threshold is determined by 
H.2(T). In addition, one notes a crossing of the curves for 
4.5 and 6 kOe, which means that H2(T) is a multivalued 
function of T . 

Figure ^| is a contour plot of the azimuthal splitting 
in a (H, T)-phase diagram. Here we have used the 3° 
contour, well above the zero-splitting error of 2°, as the 
criterion for the transition onset field. Below 10 K, H2(T) 
is nearly constant, equal to 2.25 kOe. This is somewhat 
higher than the litterature values 0, but as shown by 
Gammel et al. B, H2 depends strongly on the sample 
purity as this affects the range of the nonlocal flux line 
interactions. Above 10 K, H2{T) increases sharply, bend- 
ing away from H C 2 (T) and thereby avoiding to cross the 
upper critical field. From this it is clear that the square 
FLL is not stable at H C 2(T). The strong temperature de- 
pendence contradicts predictions based on extended GL 
theory 013), where H 2 (T) intercepts H c2 (T) at a finite 
field. An extrapolation of our measurements suggest that 
in the vicinity of H C 2(T), the FLL may become hexago- 
nal. 

Theoretically the nonlocal contribution to the vortex- 
vortex interaction weakens when the temperature rises, 
and H2(T) should somewhat increase. This relatively 
weak increase is essentially captured in the London ap- 
proach, where the order parameter |A| is taken as a con- 
stant in space Km. Moreover, |A| is assumed equal to 
the zero-field uniform value Ao (T) which determines the 
condensation energy — Fs of uniform superconductors 
(oc Aq at low temperatures and oc Ag near T c ). Being 
constant in space and field independent, the condensation 
energy so defined is commonly disregarded in calculations 
of the equilibrium FLL. Within this scheme, |A| affects 
only the value of the penetration depth, A cx 1/|A|. The 
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model is good for flux line spacings large relative to the 
core size i.e. in small fields, and for temperatures 

well below T c . 

However, with increasing density of vortices at a hxed 
T, the spatial average (|A|) decreases thus causing an 
overall increase of the system energy p5[ . The (|A|) 
suppression depends on the ratio a/£(T), with a being 
the intervortex spacing. Therefore, at a given flux line 
density, B/fo, the system energy can be reduced if the 
intervortex distance is maximized. This can be inter- 
preted as an extra (to the London interaction) repulsion 
of vortices. This repulsion is isotropic for cubic materials 
or for tetragonal crystals with H || c, and will lead to 
a deviation from square symmetry and eventually to a 
hexagonal FLL (or triangular with the structure deter- 
mined only by the ratio £ a /£b for H | c). We speculate 
that this effect should become dominant on approaching 
H C 2(T) where |A| is strongly suppressed, in accord with 
the data presented in this letter. 

In addition, we have earlier reported that in 
TmNi2B2C at low temperatures the FLL is hexagonal 
near H C 2 (T) . With decreasing field, the FLL transforms 
first to a rhombic structure followed by the square [ jl6| . 
As H is further reduced, the FLL should go through the 
same evolution in inverse order, confirmed by low field 
decoration data p| . This material is antiferromagnetic 
below 1.5 K, and it is possible that the magnetic order- 
ing has a detrimental effect on the nonlocality range in 
fields near H C 2(T). If so, both the antiferromagnetism 
and the order parameter suppression would have similar 
effects on the FLL, which may explain why the domain 
of "reverse evolution" is so broad in TmNi2B2C. Finally, 
also in agreement with our hypothesis are data on the 
FLL in VsSi that show how a distorted triangular lattice 
which is seen at T < 5 K and H = 10 kOe evolves toward 
hexagonal when T approaches T C (H) |T?|) . 

Within the microscopic theory applied to the mixed 
state, there is no artificial separation of the electromag- 
netic (London) intervortex forces and those due to the 
varying order parameter. This problem, however, has 
never been addressed for the FLL |r| . The GL approach 
does provide such a complete description, but it is valid 
only near the zero-field T c q. The accuracy of results ob- 
tained within GL for T < T c q is difficult to control. More- 
over, the very existence of the GL expansion in powers 
of the order parameter gradients away from T c q is ques- 
tionable. The coherence length £(T) does not diverge as 
T — > T c (H) ; as a consequence the gradients of A are not 
necessarily small, and the convergence of the GL series 
is difficult (if not impossible) to monitor. In this respect, 
calculations based upon keeping certain quartic terms in 

can only 



hexagonal symmetry transition in LuN^L^C. Our data 
show that the transition field H2(T), bends away from 
the upper critical field line, iJ c2 (T). This suggests that 
while approaching the H C 2(T) curve, the FLL is driven 
towards hexagonal symmetry. We propose a possible 
qualitative argument for why this might be the case. This 
poses new questions and challenges concerning the phase 
diagram of the flux line solid. 
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be justified by a microscopic theory, which is still to be 
done. 

In conclusion, we have in this letter reported measure- 
ment of the temperature dependence of the FLL square to 
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FIG. 1. FLL diffraction patterns from LuNi2B2C obtained 
at temperatures of 10 and 12 K and fields of 1.5 and 3 kOe. 
The orientation relative to the crystallographic axes, and the 
FLL (l,0)-peak azimuthal width is shown in the top panel. 
Each diffraction pattern is a sum of 4 measurements with the 
sample oriented to satisfy the Bragg condition for each of the 
main peaks/groups of peaks, and subtracted by background 
measurements taken at 20 K. The neutron wavelengths used 
were respectively A n = 15.6 A (1.5 kOe) and 11.0 A (3 kOe). 
The scattering close to q — is an artifact of imperfect back- 
ground subtraction. 

FIG. 2. Top panel: The FLL (l,0)-reflection azimuthal 
width obtained by fitting a single Gaussian to the azimuthal 
intensity distribution, as a function of temperature at an ap- 
plied field of 3 kOe. Inset: FLL azimuthal intensity distribu- 
tion folded into one quadrant, at 3 kOe and 2 and 12 K. The 
2 K data (offset for clarity) are well fitted by a single Gaus- 
sian of width 13.8° equal to the experimental resolution. The 
12 K data are fitted by two equal weight Gaussians with the 
same width, yielding a splitting of 10.5° ±0.4°. Bottom panel: 
Azimuthal splitting versus temperature for applied fields be- 
tween 1.5 and 6 kOe. The error bars includes both the error 
on the direct fit of the splitting, and the uncertainty in the 
experimental resolution. The dashed lines are guides to the 
eye. 

FIG. 3. (H, T)-phase diagram showing constant azimuthal 
splitting contours. The dashed line corresponds to a splitting 
of 3° , and is taken as the FLL symmetry transition onset field 
H2(T). Approaching H C 2(T) the scattered intensity vanishes 
and the shaded area shows the range of our measurements. A 
schematic illustration of the FLL symmetry and apex angle 
is shown for the square (centre) and hexagonal case (low H, 
high T). 
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